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Abstract 

The process pp -> W ± Z + X ifv^t^ + X is calculated to 0(a s ) 
for general C and P conserving WW Z couplings. At the Tevatron center of 
mass energy, the QCD corrections to WZ production are modest. At Large 
Hadron Collider (LHC) energies, the inclusive QCD corrections are large, but 
can be reduced significantly if a jet veto is imposed. Sensitivity limits for 
the anomalous WW Z couplings are derived from the next-to-leading order 
Z boson transverse momentum distribution for Tevatron and LHC energies. 
Unless a jet veto is imposed, next-to-leading order QCD corrections decrease 
the sensitivity to anomalous WW Z couplings considerably at LHC energies, 
but have little influence at the Tevatron. We also study, at next-to-leading 
order, rapidity correlations between the W and Z decay products, and the 
ZZ/W Z and WZ/Wj cross section ratios. These quantities are found to be 
useful tools in searching for the approximate zero present in the Standard 
Model WZ helicity amplitudes. The prospects for observing the approximate 
amplitude zero at the Tevatron and the LHC are critically assessed. 
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I. INTRODUCTION 



The electroweak Standard Model (SM) based on an SU(2) ® U(l) gauge theory has been 
remarkably successful in describing contemporary high energy physics experiments. The 
three vector boson couplings predicted by this non-abelian gauge theory, however, remain 
largely untested. The production of WZ pairs at hadron colliders provides an excellent 
opportunity to study the WW Z vertex |IS|^,[|,|5J . In addition, the reaction pp — > W ± Z 
is of interest due to the presence of an approximate zero in the amplitude of the parton 
level subprocess qiq2 — > W^Z || in the SM, which is similar in nature to the well-known 
radiation zero in the reaction pp — > W ±, y ||. In the SM, the WW Z vertex is completely 
fixed by the SU(2) ®U(1) gauge structure of the electroweak sector. A measurement of the 
WW Z vertex thus provides a stringent test of the SM. 

WZ production at hadron colliders has recently received attention due to the observation 
of a clean W + Z — > e + u e e + e~ candidate event by the CDF Collaboration 0. Double leptonic 
WZ decays are relatively background free and therefore provide an excellent testing ground 
for anomalous WW Z couplings. With an integrated luminosity of the order of 1 fb -1 , which 
is envisioned for the Main Injector era ||, a sufficient number of events should be available 
to commence a detailed investigation of the WW Z vertex in the W ± Z — > iffii^i^ channel 
(£i, £ 2 = e, /i). The prospects for a precise measurement of the WW Z couplings in this 
channel would further improve if integrated luminosities on the order of 10 fb -1 could be 
achieved (a luminosity upgraded Tevatron will henceforth be denoted by TeV*) and/or the 
energy of the Tevatron could be doubled to \fs = 3.5 — 4 TeV (an energy upgraded Tevatron 
will henceforth be referred to as the DiTevatron) ||. At the Large Hadron Collider (LHC, 
pp collisions at \/s = 14 TeV 0), it should be possible to determine the WW Z couplings 
with high precision 0. 

In contrast to low energy data and high precision measurements at the Z peak, collider 
experiments offer the possibility of a direct, and essentially model independent, determina- 
tion of the three vector boson vertices. Hadronic production of WZ pairs was first calculated 
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in Ref. M. The 0(a s ) QCD corrections to the reaction pp — ► W ± Z were first evaluated 
in Refs. JTO] and [O. Studies on the potential for probing the WWZ vertex have been 
performed in Refs. |§ and Q. 

Previous studies on probing the WWZ vertex via hadronic WZ production have been 
based on leading-order (LO) calculations HU]. In general, the inclusion of anomalous cou- 
plings at the WW Z vertex yields enhancements in the WZ cross section, especially at large 
values of the W or Z boson transverse momentum, p T (W) or p T (Z), and at large values of 
the WZ invariant mass, Mwz- Next-to-leading-order (NLO) calculations of hadronic WZ 
production have shown that the 0(a s ) corrections are large in precisely these same regions 
[|T0| , [TT|1 . It is thus vital to include the NLO corrections when using hadronic WZ production 
to test the WWZ vertex for anomalous couplings. 

In this paper, we calculate hadronic WZ production to 0(a s ), including the most general, 
C and P conserving, anomalous WWZ couplings. Our calculation also includes the leptonic 
decays of the W and Z bosons in the narrow width approximation. Decay channels where 
the W or Z boson decays hadronically are not considered here. The calculation has been 
performed using the Monte Carlo method for NLO calculations ||12|| . With this method, it 
is easy to calculate a variety of observables simultaneously and to implement experimental 
acceptance cuts in the calculation. It is also possible to compute the NLO QCD corrections 
for exclusive channels, e.g., pp> — > WZ + jet. Apart from anomalous contributions to the 
WWZ vertex we assume the SM to be valid in our calculation. In particular, we assume the 
couplings of the W and Z bosons to quarks and leptons are as given by the SM. Section II 
briefly summarizes the technical details of our calculation. 

The results of our numerical simulations are presented in Sees. Ill and IV. In contrast 
to the SM contributions to the gi<j2 — > WZ helicity amplitudes, terms associated with non- 
standard WW Z couplings grow with energy. The WZ invariant mass distribution, the 
cluster transverse mass distribution, and the Z boson transverse momentum spectrum are 
therefore very sensitive to anomalous WWZ couplings. In Sec. Ill, we focus on the impact 
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QCD corrections have on these distributions, and the sensitivity limits for the anomalous 
WW Z couplings which can be achieved at the Tevatron, DiTevatron, and LHC with various 
integrated luminosities. At LHC energies, the inclusive NLO QCD corrections in the SM 
are found to be very large at high p T (W) or p T (Z), and have a non-negligible influence on 
the sensitivity bounds which can be achieved for anomalous WW Z couplings. The large 
QCD corrections are caused by the combined effects of destructive interference in the Born 
subprocess, a log-squared enhancement factor in the q±g — > WZqi partonic cross section 



at high transverse momentum llljl , and the large quark-gluon luminosity at supercollider 
energies. At the Tevatron, on the other hand, the 0(a s ) QCD corrections are modest and 
sensitivities are only slightly affected by the QCD corrections. In Sec. Ill, we also show 
that the size of the QCD corrections at high p T (Z) or p T (W) can be significantly reduced, 
and a significant fraction of the sensitivity lost at LHC energies can be regained, if a jet 
veto is imposed, i.e., if the WZ + jet exclusive channel is considered. We also find that 
the residual dependence of the NLO WZ + jet cross section on the factorization scale 
Q 2 is significantly smaller than that of the 0(a s ) cross section for the inclusive reaction 
pp ^WZ + X. 

In Sec. IV, we study rapidity correlations between the W and Z decay products, the 
ZZ/WZ and WZ/Wj cross section ratios, and how these quantities are affected by QCD 
corrections. In the SM at tree level, the process qiq2 — > WZ exhibits an approximate 
amplitude zero at cos 6* « ^tan 2 ^ ~ 0.1 (cos 6* « — ^tan 2 ^ ~ —0.1) for ud — > 
W + Z (du — > W~Z) || which is similar in nature to the well-known radiation zero in W^ 
production in hadronic collisions. Here G* is the Z boson scattering angle in the parton 
center of mass frame relative to the quark direction and #w is the weak mixing angle. 
The radiation zero in Wj production can be observed rather easily at the Tevatron in the 
distribution of the rapidity difference Ay(7,f) = 7/(7) — y(£) of the photon and the charged 



lepton which originates from the W decay, W — > Iv |Tj|. In the SM, the Ay (7, €) distribution 
exhibits a pronounced dip at Ay(j,£) « =p0.3 in W ±r j production, which originates from 
the radiation zero. We find that, at leading order, the approximate amplitude zero in 



pp — > W ± Z — > ifi/ii^i^ leads to a dip in the corresponding Ay(Z,£i) = y(Z) — y(£i) 
distribution which is located at Ay(Z,£i) « ±0.5 (= 0) in (pp) collisions. NLO QCD 
corrections tend to fill in the dip; at LHC energies they obscure the signal of the approximate 
amplitude zero almost completely, unless a jet veto is imposed. 

Alternatively, cross section ratios can be considered. Analogous to the Z7/W7 cross 
section ratio |[L4j| , the ZZ/WZ cross section ratio as a function of the minimum transverse 
momentum of the Z boson is found to reflect the approximate amplitude zero. The ratio 
of the WZ to Wj cross sections, on the other hand, measures the relative strength of the 
radiation zero in W~f production and the approximate zero in qiq 2 — > WZ. QCD corrections 
significantly affect the ratio of ZZ to WZ cross sections, but largely cancel in the WZ/W'y 
cross section ratio. Although rapidity correlations and cross section ratios are useful tools 
in searching for the approximate amplitude zero in WZ production, it will not be easy to 
establish the effect at the Tevatron or LHC, due to the limited number of W ± Z — > if vxl^l^ 
events expected. Our conclusions, finally, are given in Sec. V. 



II. CALCULATIONAL TOOLS 



Our calculation generalizes the results of Ref. |15[ to arbitrary C and P conserving WW Z 
couplings. It is carried out using a combination of analytic and Monte Carlo integration 
techniques. Details of the method can be found in Ref. [0. The calculation is done using 



the narrow width approximation for the leptonically decaying W and Z bosons. In this 
approximation it is particularly easy to extend the NLO calculation of WZ production with 
on-shell W and Z bosons to include the leptonic W and Z decays. Furthermore, non- 
resonant Feynman diagrams such as du — ► W~* — ► e~v e Z followed by Z — > fi + contribute 
negligibly in this limit, and thus can be ignored. 
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A. Summary of 0{a s ) WZ Production including W and Z decays 



The NLO calculation of WZ production includes contributions from the square of the 
Born graphs, the interference between the Born graphs and the virtual one-loop diagrams, 
and the square of the real emission graphs. The basic idea of the method employed here is to 
isolate the soft and collinear singularities associated with the real emission subprocesses by 
partitioning phase space into soft, collinear, and finite regions. This is done by introducing 
theoretical soft and collinear cutoff parameters, 5 S and 5 C . Using dimensional regularization, 
the soft and collinear singularities are exposed as poles in e (the number of space-time 
dimensions is N = 4 — 2e with e a small number). The infrared singularities from the 
soft and virtual contributions are then explicitly canceled while the collinear singularities 
are factorized and absorbed into the definition of the parton distribution functions. The 
remaining contributions are finite and can be evaluated in four dimensions. The Monte Carlo 
program thus generates n-body (for the Born and virtual contributions) and (n + l)-body 
(for the real emission contributions) final state events. The n- and (n + l)-body contributions 
both depend on the cutoff parameters S s and 8 C , however, when these contributions are added 
together to form a suitably inclusive observable, all dependence on the cutoff parameters 
cancels. The numerical results presented in this paper are insensitive to variations of the 
cutoff parameters. 

Except for the virtual contribution, the 0(a s ) corrections are all proportional to the Born 
cross section. It is easy to incorporate the decays W — > l\V\ and Z — > £2^2 m to those terms 
which are proportional to the Born cross section; one simply replaces da Born (qiq2 — * WZ) 
with da Born (qiq2 — > WZ —>■ txVxi'^t'^) in the relevant formulae. When working at the ampli- 
tude level, the W and Z decays are trivial to implement; the vector boson polarization vector 
e ^(^)> V = W, Z , is simply replaced by the respective decay current JY(k) in the amplitude. 
Details of the amplitude level calculations for the Born and real emission subprocesses can 
be found in Ref. filf] . For £\ = £ 2 the amplitudes in principle should be antisymmetrized. 
Since the leptons originating from the decay of the W and Z bosons are usually well sepa- 
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rated, effects from the antisymmetrization of the amplitudes are expected to be very small 
and hence are ignored here. 

The only term in which it is more difficult to incorporate the W and Z decays is the 
virtual contribution. Rather than undertake the non-trivial task of recalculating the virtual 
correction term for the case of leptonically decaying weak bosons, we have instead opted to 
use the virtual correction for real on-shell W and Z bosons which we subsequently decay 
ignoring spin correlations. When spin correlations are ignored, the spin summed squared 
matrix element factorizes into separate production and decay squared matrix elements. 

Neglecting spin correlations slightly modifies the shapes of the angular distributions of 
the final state leptons, but does not alter the total cross section as long as no angular cuts 
(e.g., rapidity cuts) are imposed on the final state leptons. For realistic rapidity cuts, cross 
sections are changed by typically 10% if spin correlations are neglected. Since the size of 
the finite virtual correction is less than ~ 10% the size of the Born cross section, the overall 
effect of neglecting the spin correlations in the finite virtual correction is expected to be 
negligible compared to the combined 10% ~ 20% uncertainty from the parton distribution 
functions, the choice of the scale Q 2 , and higher order QCD corrections. 



B. Incorporation of Anomalous WW Z Couplings 



The WWZ vertex is uniquely determined in the SM by SU(2)®U(l) gauge invariance. 
In WZ production the W and Z bosons couple to essentially massless fermions, which 
insures that effectively d^W = and d^Z^ = 0. These conditions together with Lorentz 
invariance and conservation of C and P, allow three free parameters, g%,K, and A, in the 
WWZ vertex. The most general Lorentz, C, and P invariant vertex is described by the 



effective Lagrangian [17 



WWZ 



-i e cot9 



w 



giiw^w^z" - wlz v w^) + kwIw v z^ 



-^-w\ W^Z uX 
m w 



(1) 
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where and Z 11 are the W~ and Z fields, respectively, W^ v = d^W v — d v W^ and 
Z^v = d fM Z u — d v Z^. At tree level in the SM, gi — 1, k — 1, and A = 0. All higher 
dimensional operators are obtained by replacing V M with (d 2 ) m V tl , V = W, Z, where m is 
an arbitrary positive integer, in the terms proportional to Agi = g± — 1, An — k — 1, and 
A. These operators form a complete set and can be summed by replacing Agi, An, and A 
by momentum dependent form factors. All details are contained in the specific functional 
form of the form factor and its scale App. The form factor nature of Agi, An, and A will 
be discussed in more detail later in this section. 



Following the standard notation of Ref. ||17|| , we have chosen, without loss of generality, 



the W boson mass, M\y, as the energy scale in the denominator of the term proportional to 
A in Eq. ([!]). If a different mass scale, M, had been used, then all of our subsequent results 
could be obtained by scaling A by a factor M 2 /M^. 

At present, the WW Z coupling constants are only weakly constrained experimentally. 
The CDF Collaboration recently presented preliminary 95% confidence level (CL) limits on 
Ak and A from a search performed in the pp — > WZ — > jj£ + £~ (£ = e, /i) channel at large 
di-jet transverse momenta, Pr(jj) > 100 GeV [0: 

-8.6 <Ak< 9.0 (for A = 0) , -1.7 < A < 1.7 (for Ak = 0) . (2) 

Assuming SM WW'j couplings, CDF also obtained a limit on An from the reactions pp — ► 
W + W~, W ± Z -> e ± ujj with p T (jj) > 130 GeV 00: 



-1.3 < Ak < 1.4 (for A = 0) . (3) 

To derive these limits, a dipole form factor with scale A^ = 1.5 TeV was assumed (see 
below), however, the experimental bounds are quite insensitive to the value of App. Although 
bounds on these couplings can also be extracted from low energy data and high precision 
measurements at the Z pole, there are ambiguities and model dependencies in the results 
n|20||n| |. From loop contributions to rare meson decays such as Kl — > fi + fj>~ \T1 



or 



/ i+ A t [Hi! e V e HH j an d the Z — > bb width ||25|| , one estimates limits for the 
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non-standard WWZ couplings of (9(1 — 10). No rigorous bounds can be obtained from 
LEP I data if correlations between different contributions to the anomalous couplings are 
taken fully into account. In contrast, invoking a "naturalness" argument based on chiral 
perturbation theory P5,27|, one expects deviations from the SM of O(10~ 2 ) or less for gi, 



k, and A. 

If C or P violating WWZ couplings are allowed, four additional free parameters, g^, 
#5, R, and A appear in the effective Lagrangian. For simplicity, these couplings are not 
considered in this paper. 

The Feynman rule for the WWZ vertex factor corresponding to the Lagrangian in Eq. ([[]) 

is 



i 9wwz Qw r^jX&j ki, k 2 ) — —i gwwz Qw 



rj^(Mi,*2) + r^(Mi,k) 



(4) 



where the labeling conventions for the four-momenta and Lorentz indices are defined by 
Fig. [1|, gwwz = e cot # w is the WWZ coupling strength, Q w is the electric charge of the 
W boson in units of the proton charge e, and the factors T SM and r NSM are the SM and 
non-standard model vertex factors: 



T f^v( k i k i,k 2 ) = {h- k 2 ) p g Ufl + 2k fl gp v -2k v g Pfl , (5) 
Tf^(k, h, k 2 ) = 1 U 9l + + A J^) (h - k 2 )f3 g Vi x (6) 



w , 



(k 1 -k 2 )pk I/ k fl + (Ag 1 + AK + X)k^ l gp v 



M w 



The non-standard model vertex factor is written here in terms of Agi = g\ — 1, Ak = k — 1, 
and A, which all vanish in the SM. 

It is straightforward to include the non-standard model couplings in the amplitude level 
calculations. Using the computer algebra program FORM | |2^| , we have computed the 



qiq 2 — > WZ virtual correction with the modified vertex factor of Eq. (HI), however, the 
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resulting expression is too lengthy to present here. The non-standard WW Z couplings of 
Eq. ([1]) do not destroy the renormalizability of QCD. Thus the infrared singularities from 
the soft and virtual contributions are explicitly canceled, and the collinear singularities are 
factorized and absorbed into the definition of the parton distribution functions, exactly as 
in the SM case. 

The anomalous couplings can not be simply inserted into the vertex factor as constants 
because this would violate S'-matrix unitarity. Tree level unitarity uniquely restricts the 



WW Z couplings to their SM gauge theory values at asymptotically high energies |29j . This 
implies that any deviation of Api, Ak, or A from the SM expectation has to be described by a 
form factor Ag^M^ z , p 2 w , p 2 z ) , Ak(M^ z ,p^,p 2 z ), or A(M^ z , p 2 w , p|) which vanishes when 
either the square of the W Z invariant mass, M^ z , or the square of the four-momentum 
of the final state W or Z {ply or p 2 z ) becomes large. In WZ production ply m Myy and 
p\ ~ M§ even when the finite W and Z widths are taken into account. However, large values 
of Mfyz wm De probed at future hadron colliders like the LHC and the M^ z dependence 
of the anomalous couplings has to be included in order to avoid unphysical results which 
would violate unitarity. Consequently, the anomalous couplings are introduced via form 
factors PD|PT[ 



Af7° 

&9i(M 2 wz ,p 2 w = M 2 w ,p 2 z = Ml) = n ../i , (7) 

[i -\- ivi wz / i\ FF ) 

Ak{M 2 wz , P 2 w = M 2 w ,p 2 z = Ml) = n . A ^ K ° /A2 , n , (8) 

I J- "I" m wzl iv ff) 

X(M^ z ,p 2 w = M 2 w ,p\ = Mf) = 7——^ /A2 , n , (9) 

I I + m WZl ly FF) 

where Ag®, Ak°, and A are the form factor values at low energies and Aff represents the 
scale at which new physics becomes important in the weak boson sector, e.g., due to new 
resonances or composite structures of the W and Z bosons. In order to guarantee unitarity, 
it is necessary to have n > 1/2 for An and n > 1 for Ag± and A. For the numerical 
results presented here, we use a dipole form factor {n = 2) with a scale Aff = 1 TeV, 
unless explicitly stated otherwise. The exponent n = 2 is chosen in order to suppress WZ 
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production at energies \/I > Aff ^> Mw, Mz, where novel phenomena like resonance or 
multiple weak boson production are expected to become important. 

Form factors are usually not introduced if an ansatz based on chiral perturbation theory 
is used. In the framework of chiral perturbation theory, the effective Lagrangian describing 
the anomalous vector boson self-interactions breaks down at center of mass energies above 
a few TeV p^ , p7[] (typically A-nv ~ 3 TeV, where v ~ 246 GeV is the Higgs field vacuum 
expectation value). Consequently, one has to limit the center of mass energies to values 
sufficiently below 4ttv in this approach. 

III. QCD CORRECTIONS AND NON-STANDARD WW Z COUPLINGS 

We shall now discuss the phenomenological implications of NLO QCD corrections to W Z 
production at the Tevatron (pp collisions at y/s = 1.8 TeV) and the LHC (pp collisions at 

= 14 TeV). We also consider a possible upgrade || of the Tevatron to = 3.5 TeV 
(DiTevatron). We first briefly describe the input parameters, cuts, and the finite energy 
resolution smearing used to simulate detector response. We then discuss in detail the impact 
of NLO QCD corrections on the observability of non-standard WW Z couplings in W Z 
production at the Tevatron, DiTevatron, and LHC. To simplify the discussion, we shall 
concentrate on W + Z production. In pp collisions the rates for W + Z and W~~ Z production 
are equal. At pp colliders, the W~Z cross section is about 30% smaller than that of W + Z 
production. Furthermore, we shall only consider W + —>■ i\v\ and Z — > decays (£i, £2 = 
e, /i). 

A. Input Parameters 

The numerical results presented here were obtained using the two-loop expression for 
a s . The QCD scale Aqcd is specified for four flavors of quarks by the choice of the parton 
distribution functions and is adjusted whenever a heavy quark threshold is crossed so that 
a s is a continuous function of Q 2 . The heavy quark masses were taken to be m?, = 5 GeV 
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and m t = 150 GeV, which is consistent with the bound obtained by D0, mt > 131 GeV [32 
and the value suggested by the current CDF data, m t = 174 ± 10 IjJ GeV Our results 
are insensitive to the value chosen for mt- 

The SM parameters used in the numerical simulations are Mz = 91.173 GeV, My/ = 
80.22 GeV, a{M w ) = 1/128, and sin 2 6 W = 1-(M W /M Z ) 2 . These values are consistent with 



recent measurements at LEP, SLC, the CERN pp collider, and the Tevatron p4] , |3^ . |36| . |37 . 
The soft and collinear cutoff parameters, as discussed in Sec. IIA, are fixed to 5 S = 10~ 2 and 
5 C = 10~ 3 . The parton subprocesses have been summed over u, d, s, and c quarks and the 
Cabibbo mixing angle has been chosen such that cos 2 9c = 0.95. The leptonic branching 
ratios have been taken to be B(W -> Iv) = 0.107 and B(Z -> g+f) = 0.034 and the total 
widths of the W and Z bosons are Ty^ = 2.12 GeV and Y z = 2.487 GeV. Except where 
otherwise stated, a single scale Q 2 = M^ z , where M W z is the invariant mass of the WZ 
pair, has been used for the renormalization scale /x 2 and the factorization scale M 2 . 

In order to get consistent NLO results it is necessary to use parton distribution functions 
which have been fit to next-to-leading order. In our numerical simulations we have used 
the Martin-Roberts-Stirling (MRS) |3§ set SO' distributions with A 4 = 230 MeV, which 



take into account the most recent NMC [39| and CCFR [[5J data and are consistent with 



measurements of the proton structure functions at HERA ||41|| . For convenience, the MRS 



set SO' distributions have also been used for the LO calculations. 

B. Cuts 

The cuts imposed in our numerical simulations are motivated by the finite acceptance of 
detectors. The complete set of transverse momentum (p T ), pseudorapidity (77), and separa- 
tion cuts can be summarized as follows. 
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Tevatron 



LHC 



p T (£) > 20 GeV p T (£) > 25 GeV 

i> T > 20 GeV i> T > 50 GeV 

< 2.5 \r}(£)\ < 3.0 

AR(£,£)>0A AR(£,£)>0A 

Here, AR = [(A0) 2 + (A?]) 2 ] 1 / 2 is the separation in the pseudorapidity - azimuthal angle 
plane. The AR(£, £) cut is only imposed on leptons of equal electric charge. It has only a 
small effect on the WZ cross section. For simplicity, identical cuts are imposed on final state 
electrons and muons. The large missing transverse momentum {j) T ) cut at LHC energies, 
which severely reduces the total WZ cross section, has been chosen to reduce potentially 
dangerous backgrounds from event pileup J42|, pp — > Zbb — > ^vi^^ + X, and processes 



where particles outside the rapidity range covered by the detector contribute to the missing 



transverse momentum. Present studies f43| , f44|1 indicate that these backgrounds are under 



control for j> T > 50 GeV. The total W + Z cross section within cuts in the Born approximation 
at the Tevatron, DiTevatron, and LHC is 8.5 fb, 22.4 fb, and 25.9 fb, respectively. If the 
LHC is operated significantly below the design luminosity of L = 10 34 cm~ 2 s _1 0, the 
background from event pileup is less severe and it may well be possible to lower the missing 
p T threshold. If the j> T > 50 GeV cut is replaced by a j> T > 20 GeV requirement, the total 
LO W + Z cross section triples to 80.8 fb. 



C. Finite Energy Resolution Effects 

Uncertainties in the energy measurements of the charged leptons in the detector are 
simulated in our calculation by Gaussian smearing of the particle four-momentum vector 
with standard deviation a. For distributions which require a jet definition, e.g., the WZ + 
1 jet exclusive cross section, the jet four-momentum vector is also smeared. The standard 
deviation a depends on the particle type and the detector. The numerical results presented 
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here for the Tevatron/DiTevatron and LHC center of mass energies were made using a values 
based on the CDF ||5| and ATLAS [j43] specifications, respectively. 



D. Inclusive NLO Cross Sections 

The sensitivity of WZ production to anomalous WW Z couplings in the Born approxi- 
mation was studied in detail in Refs. and Q. The distributions of the Z boson transverse 
momentum, p T (Z), and the WZ invariant mass, M WZl were found to be sensitive to the 
anomalous couplings. However, at hadron colliders the WZ invariant mass cannot be de- 
termined unambiguously because the neutrino from the W decay is not observed. If the 
transverse momentum of the neutrino is identified with the missing transverse momentum 
of a given WZ event, the unobserved longitudinal neutrino momentum Pl{v) can be recon- 
structed, albeit with a twofold ambiguity, by imposing the constraint that the neutrino and 
the charged lepton four-momenta combine to form the W rest mass [ IB, 17 . Neglecting the 
lepton mass one finds 

Pl{u) = 2pW) ( PlW ( M ^ + 2Pt{£) ' *r) (10) 



±p(£) 



(M w + 2p T (£)-^ T ) 2 -Ap 2 T {£)i>, 



nl/2- 



where Pl(£) denotes the longitudinal momentum of the charged lepton. The two solutions 
for p L {y) are used to reconstruct two values for M wz . Both values are then histogrammed, 
each with half the event weight. 

The differential cross section for Mwz in the reaction pp — > W + Z+X — > l\V\l\l2 +A at 
^/s = 1.8 TeV is shown in Fig. |^. The Born and NLO results are shown in Fig. ^|a and Fig. 0b, 
respectively. In both cases, results are displayed for the SM and for three sets of anomalous 
couplings, namely, (A = -0.5, Ak° = Ag® = 0), (Ak° = -1.0, A = A#° = 0), and (A#° = 
—0.5, Aft = A = 0). For simplicity, only one anomalous coupling at a time is allowed 
to differ from its SM value. The figure shows that at the Tevatron center of mass energy, 
NLO QCD corrections do not have a large influence on the sensitivity of the reconstructed 
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WZ invariant mass distribution to anomalous couplings. The 0(a s ) corrections at Tevatron 
energies are approximately 30% for the SM as well as for the anomalous coupling cases. Since 
the anomalous terms in the helicity amplitudes grow like y§/Mw (s/Myy) for An (A and 
Agi) ||, where s denotes the parton center of mass energy squared, non-standard couplings 
give large enhancements in the cross section at large values of Mwz- 

The WZ invariant mass distributions at the DiTevatron and the LHC are shown in Figs. ||| 
and U, respectively. In both cases, the sensitivity of the Mwz distribution to anomalous 
WWZ couplings is significantly more pronounced than for pp collisions at = 1.8 TeV. 
NLO QCD corrections enhance the SM Mwz differential cross section by about a factor 2 
at the LHC, whereas the 0(a s ) corrections at the DiTevatron are very similar in size to 
those found at Tevatron energies. For non-standard WWZ couplings, the QCD corrections 
are more modest at the LHC. Because of the form factor parameters assumed, the result 
for Ak° = —1 approaches the SM result at large values of Mwz- As mentioned before, we 
have used n = 2 and a form factor scale of App = 1 TeV in all our numerical simulations 
[see Eqs. (|7|) - (0)]. For a larger scale Aff, the deviations from the SM result become more 
pronounced at high energies. No significant change in the shape of the Mwz distribution is 
observed for all center of mass energies considered. 

In Fig. H we investigate the influence of anomalous WWZ couplings on the WZ invariant 
mass spectrum at next-to-leading order, together with the effect of the twofold ambiguity in 
Mwz originating from the reconstruction of the longitudinal momentum of the neutrino from 
the W decay, in more detail. The lower dotted, dashed, and dot-dashed lines display the WZ 
invariant mass distribution for Ak° = + 1 (+1), Ag\ = +0.5 (+0.25), and A = +0.5 (+0.25), 
whereas the upper curves show the Mwz spectrum for Ak° = —1 (—1), Ag® = —0.5 (—0.25), 
and A = —0.5 (—0.25) at the Tevatron (LHC). For Ak° and Ag®, negative anomalous 
couplings lead to significantly larger deviations from the SM prediction than positive non- 
standard couplings of equal magnitude, whereas there is little difference for A (dashed lines). 
For Ak° the sign dependence is more pronounced at small energies. Other distributions, 
such as the cluster transverse mass distribution or the transverse momentum distribution of 
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the Z boson display a similar behaviour. 

This effect can be easily understood from the high energy behaviour of the WZ produc- 
tion amplitudes, A4(Xz, Xw), where Xz (Xw) denotes the helicity of the Z (W) boson HQ. 
In the SM, only .M(±, =p) and Ai(0, 0) remain finite for s — > oo. Contributions to the he- 
licity amplitudes proportional to A mostly influence the (±, ±) amplitudes and increase like 
s/M w at large energies. The SM (±, ±) amplitudes vanish like 1/s, and the non-standard 
terms dominate except for the threshold region, ~ Mw + Mz- For non-standard values of 
A, the cross section therefore depends only very little on the sign of the anomalous coupling. 
Terms proportional to Agi also increase like s/M w with energy, but mostly contribute to 
the (0,0) amplitude, which remains finite in the SM in the high energy limit. Interference 
effects between the SM and the anomalous contributions to the (0,0) amplitude, thus, are 
non- negligible, resulting in a significant dependence of the cross section on the sign of Agi. 
Terms proportional to An, finally, are proportional to yk/Mw and mostly influence the 
amplitudes with a longitudinal W and a transverse Z boson. In the SM, these terms vanish 
like l/vl. In the high energy limit one therefore expects little dependence of the cross 
section on the sign of Ak, similar to the situation encountered for A (see Fig. ^o). However, 
since the terms proportional to An increase less drastically with energy, interference effects 
between those terms and the SM amplitudes are substantial near threshold. 

The dash-double-dotted line in Fig. § shows the true M wz distribution. The distribution 
of the reconstructed invariant mass at both center of mass energies is harder than the true 
Mwz distribution. At LHC energies, the twofold ambiguity in the reconstructed WZ invari- 
ant mass only slightly affects the Mwz distribution (Fig. |5|b). At the Tevatron, on the other 
hand, the true and reconstructed invariant masses are quite different for Mwz > 500 GeV, 
thus severely degrading the sensitivity to non-standard WWZ couplings. If the W decay is 
treated in the narrow width approximation, one of the two reconstructed invariant masses 
coincides with the true WZ invariant mass. Since the Mwz spectrum is steeply falling, 
the incorrect solution of the reconstructed invariant mass influences the distribution in a 
noticeable way only if it is larger than the true WZ invariant mass. The average difference 
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(absolute value) between the two reconstructed values of Mwz is almost independent of the 
center of mass energy. As a result, the twofold ambiguity in the reconstructed WZ invariant 
mass affects the Mwz spectrum at the LHC much less than at the Tevatron. 

As an alternative to the WZ invariant mass spectrum, the differential cross section of 
the cluster transverse mass M r (£ 1 £^£^;^ T ) can be studied. The cluster transverse mass 
is defined by 

2 

- |p T (c) +pi T | 2 , (11) 

where M c denotes the invariant mass of the cluster c = iil^l^- The Mt distribution at the 
Tevatron and the LHC is shown in Fig. || Since QCD corrections change its shape only 
slightly, we only show the NLO Mt distribution. At the Tevatron, the cluster transverse 
mass distribution is seen to be significantly more sensitive to anomalous couplings than the 
reconstructed WZ invariant mass distribution, in particular for A and Agi. The cluster 
transverse mass distribution for pp collisions at y/s = 3.5 TeV is qualitatively very similar 
to that obtained at Tevatron energies and is therefore not shown. 

In Figs. |7] - |9| we show the differential cross section for the transverse momentum of 
the Z boson, p T (Z). The p T (Z) spectrum is seen to be considerably more sensitive to 
non-standard WW Z couplings than the cluster transverse mass or the WZ invariant mass 
distributions. At high transverse momentum, a large enhancement of the cross section is 
observed. On the other hand, at Tevatron and DiTevatron energies, the Pt(Z) differen- 
tial cross section is smaller than predicted in the SM for Pt{Z) < 30 GeV, if anomalous 
WW Z couplings are present. Due to the relatively large j> T cut imposed, the Z boson 
transverse momentum distribution at the LHC in the Born approximation displays a pro- 
nounced shoulder at p T (Z) w 65 GeV (see Fig. |9|a). The j> T cut mostly affects the small 
Pt(Z) region (p T (Z) < 100 GeV) and therefore does not significantly reduce the sensitivity 
to non-standard WW Z couplings. Once NLO QCD corrections are taken into account this 
shoulder disappears due to contributions from the real emission subprocesses, qiq 2 — > WZg 
and q±g — > WZq 2 . 
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In contrast to the other distributions studied so far, the shape of the SM Z boson 
transverse momentum spectrum is considerably affected by NLO QCD corrections. This is 
demonstrated in detail in Fig. |X0j, where we show the ratio of the NLO and LO differential 
cross sections of the Z boson transverse momentum. At Tevatron and DiTevatron energies, 
the 0(a s ) corrections are approximately 30% for the SM, and 35 - 40% for the anomalous 
coupling cases at small p T {Z) values. In the SM case, the size of the QCD corrections in- 
creases to ~ 60% for Pt{Z) — 300 GeV at the Tevatron, and to ~ 80% at the DiTevatron. 
For non-standard couplings, on the other hand, the QCD corrections are between 20% and 
40% over the whole Pj-(Z) range plotted. This is exemplified by the dashed and dot-dashed 
lines in Fig. |T0"|a, which show the ratio of NLO to LO cross sections for A = —0.5 at the 



Tevatron and DiTevatron, respectively. At the LHC (see Fig. p70|b) , the shape of the p T (Z) 
distribution is drastically altered by the NLO QCD corrections. At p?{Z) = 800 GeV, the 
QCD corrections increase the SM cross section by about a factor 5, whereas the enhance- 
ment is only a factor 1.6 at p T (Z) = 100 GeV. In the presence of anomalous couplings, 
the higher order QCD corrections are much smaller than in the SM. In regions where the 
anomalous terms dominate, the 0(a s ) corrections are typically between 20% and 40%. This 
is illustrated by the dashed curve in Fig. |TU|b, which shows the NLO to LO cross section 
ratio for A = —0.25. At next-to-leading order, the sensitivity of the Z boson transverse 
momentum spectrum to anomalous couplings is thus considerably reduced. The transverse 
momentum distribution of the W boson in the SM exhibits a similar strong sensitivity to 
QCD corrections. Qualitatively, the change of the shape of the p?{Z) distribution in WZ 
production and of the photon transverse momentum distribution in W'y production at high 
energies [(49[ are very similar. 



E. Exclusive NLO QCD Corrections 

The large QCD corrections at high values of Pt{Z) are caused by the combined effects 
of destructive interference in the Born process, a collinear enhancement factor in the q^g — > 
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WZq 2 partonic cross section for p T (Z) S> Mw, and the large qg luminosity at LHC energies. 
In the SM, delicate cancellations between the amplitudes of the Born diagrams occur in 
the central rapidity region in WZ production. These cancellations are responsible for the 
approximate amplitude zero || and suppress the WZ differential cross section, in particular 
for large W and Z transverse momenta. In the limit p T (Z) 3> Mw, the cross section for 
qig — > \VZq2 can be obtained using the Altarelli-Parisi approximation for collinear emission. 
One finds fll| : 



da(q l9 - WZq 2 ) = da(q l9 -> q,Z) ^ In 2 (j^) , (12) 

where g w = e/ sin0\y- Thus, the quark gluon fusion process carries an enhancement factor 
ln 2 (p| 1 (Z)/M^ / ) at large values of the Z boson transverse momentum. It arises from the 
kinematic region where the Z boson is produced at large p T and recoils against the quark, 
which radiates a soft W boson collinear to the quark. Since the Feynman diagrams entering 
the derivation of Eq. flT2|) do not involve the WW Z vertex, the logarithmic enhancement 
factor only affects the SM matrix elements. At the LHC, the p T {Z) differential cross section 
obtained using Eq. flT2] ) agrees within 30% with the exact Z boson transverse momentum 
distribution for Pj>{Z) > 200 GeV [IT]. Together with the very large qg luminosity at 



supercollider energies and the suppression of the SM WZ rate at large values of the Z boson 
transverse momentum in the Born approximation, the logarithmic enhancement factor is 
responsible for the size of the inclusive NLO QCD corrections to WZ production, as well as 
for the change in the shape of the Pt(Z) distribution. The same enhancement factor also 
appears in the antiquark gluon fusion process, however, the qg luminosity is much smaller 
than the qg luminosity for large values of the Z boson transverse momentum. Since the W 
boson does not couple directly to the gluon, the process q\q 2 — > WZg is not enhanced at large 
values of the Z boson transverse momentum. Arguments similar to those presented above 
also apply to the W boson transverse momentum distribution in the limit Pj-(W) 3> Mz, 
and a relation analogous to Eq. (|I^) can be derived. 

From the picture outlined in the previous paragraph, one expects that, at next-to- leading 
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order at supercollider energies, WZ events with a high p T Z boson most of the time also 
contain a high transverse momentum jet. At the Tevatron, on the other hand, the fraction 
of high Pt(Z) WZ events with a hard jet should be considerably smaller, due to the much 
smaller qg luminosity at lower energies. The decomposition of the inclusive SM NLO Pt{Z) 
differential cross section into NLO 0-jet and LO 1-jet exclusive cross sections at the Tevatron 
and LHC is shown in Figs. [TTl a and |l2]a, respectively. The SM NLO 0-jet p T (Z) distribu- 



tion at the two center of mass energies are compared with the corresponding distributions 
obtained in the Born approximation in Figs. [TT| b and |i~2"lb. Here, a jet is defined as a quark 
or gluon with 

p T (j) > 10 GeV and \rj(j)\ < 2.5 (13) 

at the Tevatron and 

p T (j) > 50 GeV and \r]{j)\ < 3 (14) 

at the LHC. The sum of the NLO 0-jet and the LO 1-jet exclusive cross section is equal to 
the inclusive NLO cross section. The results for the NLO exclusive WZ + jet and the LO 
exclusive WZ + 1 jet differential cross sections depend explicitly on the jet definition. Only 
the inclusive NLO distributions are independent of the jet definition. It should be noted 
that the jet transverse momentum threshold can not be lowered to arbitrarily small values 
in our calculation for theoretical reasons. For transverse momenta below 5 GeV (20 GeV) 
at the Tevatron (LHC), soft gluon resummation effects are expected to significantly change 
the jet p T distribution |50| . For the jet definitions listed above (Eqs. (fC|) and (|I4"|)), these 



effects are expected to be unimportant, and therefore are ignored in our calculation. 

With the jet definition of Eq. (|T3|), the inclusive NLO cross section at the Tevatron is 
composed predominately of 0-jet events (see Fig. [TT|a). Due to the logarithmic enhancement 
factor, the 1-jet cross section becomes relatively more important at large values of the Z 
boson transverse momentum. For p T {Z) values above 200 GeV, the 1-jet cross section is 
larger than the 0-jet rate at the LHC, and dominates completely at high Pt{Z) (see Fig. |l2|a). 
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Figure [Tip compares the NLO WZ + jet cross section with the result obtained in the Born 
approximation at the Tevatron. For the jet definition chosen [see Eq. fll^DI, the results are 
almost identical over the entire transverse momentum range displayed. At the LHC, the 
NLO WZ + jet result is at most 20% smaller than the cross section obtained in the Born 
approximation (Fig. |T2]b). Note that the characteristic shoulder at p T (Z) 65 GeV in 
the Born p T {Z) distribution, which results from the large f T cut, is eliminated in the NLO 
WZ + jet differential cross section. The LO and NLO WZ + jet p T (ii^), ^y, and 2/(^1,2) 
differential cross sections also agree to better than 20% The results for DiTevatron 

energies are very similar to those obtained for the Tevatron. 

If the jet defining p T threshold is lowered to 30 GeV and jets can be identified out to 
\vU)\ = 4-5 at the LHC, the NLO WZ + jet p T {Z) differential cross section is approx- 
imately 30% smaller than the result obtained with the jet definition of Eq. ([14]). Present 
studies suggest [^3| , ^4] , pT| that jets fulfilling the criteria of Eq. (|T4] ) can be identified at the 



LHC without problems, whereas it will be difficult to reconstruct a jet with a transverse mo- 
mentum of less than 30 GeV. The pseudorapidity range covered by the LHC is not expected 
to extend beyond |?7| = 4.5. 

The results shown so far were obtained for Q 2 = M^ z . Since the WZ + 1 jet and the 
WZ + jet cross section in the Born approximation are tree level results, they are sensitive 
to the choice of the factorization scale Q 2 . Figure [13] displays the scale dependence of the 
Born, the inclusive NLO, the 0(a s ) 0-jet exclusive, and the 1-jet exclusive cross sections 
for the Tevatron and LHC center of mass energies. The total cross section for the reaction 
pp — > W + Z + X — y l\v\^l^ + X is plotted versus the scale Q. The factorization scale M 2 
and the renormalization scale \j? have both been set equal to Q 2 . 

The scale dependence of the Born cross section enters only through the Q 2 dependence 
of the parton distribution functions. The qualitative differences between the results at the 
Tevatron and the LHC are due to the differences between pp versus pp scattering and the 
ranges of the x- values probed. At the Tevatron, WZ production in pp collisions is dominated 
by valence quark interactions. The valence quark distributions decrease with Q 2 for the x- 
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values probed at the Tevatron (typically x > 0.1). On the other hand, at the LHC, sea 
quark interactions dominate in the pp process and smaller x-values are probed (typically 
x ~ 0.02). The sea quark distributions increase with Q 2 for the x- values probed at the 
LHC. Thus the Born cross section decreases with Q 2 at the Tevatron but increases with Q 2 
at the LHC. 

The scale dependence of the 1-jet exclusive cross section enters via the parton distribution 
functions and the running coupling a s (Q 2 ). Note that the 1-jet exclusive cross section is 
calculated only to lowest order and thus exhibits a strong scale dependence. The dependence 
on Q here is dominated by the scale dependence of a s (Q 2 ) which is a decreasing function 
of Q 2 . At the NLO level, the Q dependence enters not only via the parton distribution 
functions and the running coupling a s (Q 2 ), but also through explicit factorization scale 
dependence in the order a s (Q 2 ) correction terms. The NLO 0-jet exclusive cross section 
is almost independent of the scale Q. It shows a non- negligible variation with the scale 
only in the region Q < 100 GeV at the Tevatron. In the WZ + jet cross section, the 
scale dependence of the parton distribution functions is compensated by that of a s (Q 2 ) 
and the explicit factorization scale dependence in the correction terms. The Q dependence 
of the inclusive NLO cross section is dominated by the 1-jet exclusive component and is 
significantly larger than that of the NLO 0-jet cross section. 

F. Sensitivity Limits 

We now study the impact that NLO QCD corrections to WZ production have on the 
sensitivity limits for Ag®, Ak°, and A at the Tevatron, DiTevatron, and LHC. For the 
Tevatron we consider integrated luminosities of JCdt = 1 fb _1 , as envisioned for the Main 
Injector era, and 10 fb _1 (TeV*) which could be achieved through additional upgrades of 
the Tevatron accelerator complex [Q . In the case of the DiTevatron we assume an integrated 
luminosity of 10 fb _1 . For the LHC we consider integrated luminosities of JCdt = 10 fb _1 
and 100 fb _1 M. To extract limits at the Tevatron, TeV*, and DiTevatron, we shall sum 
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over both W charges. For the LHC, we only consider W + Z production. Interference effects 
between Ag°, Ak°, and A are fully incorporated in our analysis. 

To derive 95% CL limits we use the Pt{Z) distribution and perform a \ 2 test ||52|| . In 



the Born approximation, the Z boson transverse momentum distribution in general yields 
the best sensitivity bounds. Furthermore, we use the cuts summarized in Sec. IIIB, and 
the jet definitions of Eqs. fll3|) and (0). Unless explicitly stated otherwise, a dipole form 



factor (n = 2) with scale = 1 TeV is assumed. For the Tevatron with 1 fb _1 we split 
the Pt{Z) distribution into 3 bins, whereas 7 bins are used in all other cases. In each bin 
the Poisson statistics are approximated by a Gaussian distribution. In order to achieve a 
sizable counting rate in each bin, all events with p T {Z) > 60 GeV (120 GeV) at the Tevatron 
(TeV*) are collected in a single bin. Similarly, all events with p T {Z) > 180 GeV (240 GeV 
[480 GeV]) at the DiTevatron (LHC with 10 fb _1 [100 fb -1 ]) are accumulated into one bin. 
This procedure guarantees that a high statistical significance cannot arise from a single 
event at large transverse momentum, where the SM predicts, say, only 0.01 events. In order 
to derive realistic limits we allow for a normalization uncertainty of 50% in the SM cross 
section. Background contributions are expected to be small for the cuts we impose, and are 
ignored in our derivation of sensitivity bounds. 



Our results are summarized in Figs. [14] and [15], and Tables I and II. The cross section in 
each bin is a bilinear function of the anomalous couplings AkP, A , and Ag®. Studying the 
correlations in the Are - A , the An - Ag°, and the Ag® - A planes is therefore sufficient 



to fully include all interference effects between the various WW Z couplings. Figure |14 
shows 95% CL contours in the three planes for the pp collider options obtained from the 
inclusive NLO p T (Z) distribution. Table I displays the 95% CL sensitivity limits, including 
all correlations, at leading order and next-to-leading order for the three WW Z couplings 
for the process pp — > W ± Z + X -> ifv^t^ + x - At Tevatron and DiTevatron energies, 
the increase of the cross section at 0(a s ) and the reduced sensitivity at large values of the 
Z boson transverse momentum balance each other, and the limits obtained at LO and NLO 
are usually very similar. Our limits fully reflect the strong sign dependence of the differential 
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cross sections observed for Aft and Ag° (see Fig. |5|). 

With an integrated luminosity of 1 fb -1 it will not be possible to perform a very precise 
measurement of the WW Z vertex in the WZ —>■ iiUii^i^ channel at the Tevatron. For 
integrated luminosities of less than a few fb -1 , the limits which can be achieved, however, 
may be significantly improved by combining the bounds from WZ — > i^v^t^t"^ with the 
limits obtained from pp — > WZ — > &\V\jj and pp — > WZ — > jj££^2 a ^ large di-jet transverse 



momenta [p3|1 - Currently, these channels are used by the CDF Collaboration to extract 
information on the structure of the WWZ vertex [0,|l8|. Decay modes where the W or 
Z boson decays hadronically have a considerably larger branching ratio than the WZ — > 
i\V\i\i2 channel and thus yield higher rates. On the other hand, they are plagued by 
a substantial W/Z + 2 jet QCD background, which, for large integrated luminosities (> 
10 fb _1 ), will eventually limit the sensitivity of the semi-hadronic WZ decay channels to 
anomalous WWZ couplings. 

The limits which can be achieved for Ak° at the TeV* from WZ — > £x v i^2^2 are about a 
factor 1.8 better than those at the Tevatron with 1 fb _1 . The bounds on Ag° and A improve 
by a factor 2 to 2.7. Increasing the energy of the Tevatron to 3.5 TeV (DiTevatron) improves 
the limits again significantly, in particular, the bound on An . Due to the rather strong 
interference effects between the SM and the anomalous terms of the helicity amplitudes for 
Agi and A/c, the contours sometimes deviate substantially from the elliptical form naively 
expected. Furthermore, significant correlations are observed, in particular, between An 
and Ag® (see Fig. |l4|b) . The limits obtained with a 0-jet requirement imposed are virtually 



identical to those resulting from the inclusive NLO p T (Z) distribution. 



The 95% CL limit contours for the LHC are shown in Fig. |15[ Table II summarizes 
the LO and NLO sensitivity bounds which can be achieved at the LHC. At supercollider 
energies, the inclusive NLO QCD corrections in the SM are very large and drastically change 
the shape of the SM p T (Z) distribution (see Fig. |9]). As a result, NLO QCD corrections 
reduce the sensitivity to anomalous couplings by 20 - 40%. As the integrated luminosity 
increases, larger transverse momenta become accessible. The difference between the LO 
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and NLO sensitivity bounds for 100 fb 1 therefore is somewhat larger than for 10 fb 1 . 
For the parameters chosen, the inclusive NLO bounds which can be obtained from pp — > 
W + Z + X -> if + X at y/s = U TeV with JCdt = 10 fb" 1 are quite similar to 

those which are expected from the DiTevatron for W^Z production and equal integrated 
luminosity. 

As we have seen in Sec. HIE, the size of the 0(a s ) QCD corrections at the LHC can be 
significantly reduced by vetoing hard jets in the central rapidity region, i.e., by imposing 
a "zero jet" requirement and considering only the WZ + jet channel. A zero jet cut for 



example has been imposed in the CDF measurement of the ratio of W to Z cross sections [54 



and the W mass measurement |]55 |. The sensitivity limits obtained for the WZ+0 jet channel 
at NLO are 10 - 30% better than those obtained in the inclusive NLO case and are quite 
often close to those obtained from the leading order Pt{Z) distribution (see Table II and 



the dotted contours in Fig. [151). The NLO WZ + jet differential cross section is also more 
stable to variations of the factorization scale Q 2 than the Born and inclusive NLO WZ + X 
cross sections (see Fig. |13| ). The systematic errors which originate from the choice of Q 2 
will thus be smaller for bounds derived from the NLO WZ + jet differential cross section 
than those obtained from the inclusive NLO WZ + X or the Born cross section. The limits 
extracted from the WZ + jet exclusive channel depend only negligibly on the jet definition 
used. 

The bounds which can be achieved at the LHC improve by up to a factor 3 if an integrated 
luminosity of 100 fb -1 can be achieved (dot-dashed contours in Fig. Note that the Ak° 



and Ag® limits are strongly correlated in this case. This effect is due to the relatively small 
form factor scale chosen (A FF = 1 TeV), which significantly suppresses the non-standard 
terms in the helicity amplitudes at high energies. 

At Tevatron (DiTevatron) energies, the sensitivities achievable are insensitive to the exact 
form and scale of the form factor for A FF > 400 GeV (A FF > 800 GeV). At the LHC, the 
situation is somewhat different and the sensitivity bounds depend on the value chosen for 
A FF |§. This is illustrated in Table lie, where we list the bounds which can be achieved at 
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the LHC with JCdt = 100 fb _1 and a form factor scale of App — 3 TeV. The limits for the 
higher scale are a factor 1.8 to 5 better than those found for App = 1 TeV with the same 
integrated luminosity. For App > 3 TeV, the sensitivity bounds depend only marginally on 
the form factor scale, due to the very rapidly falling cross section at the LHC for parton 
center of mass energies in the multi-TeV region. The dependence of the limits on the cutoff 
scale App in the form factor can be understood easily from Fig. |9|. The improvement in 
sensitivity with increasing App is due to the additional events at large Pp{Z) which are 
suppressed by the form factor if the scale App has a smaller value. 

To a lesser degree, the bounds also depend on the power n in the form factor, which we 
have assumed to be n = 2. For example, the less drastic cutoff for n = 1 instead of n = 2 
in the form factor allows for additional high p T {Z) events and therefore leads to a slightly 
increased sensitivity to the low energy values AkP, Ag®, and A . The sensitivity bounds 
listed in Tables I and II can thus be taken as representative for a wide class of form factors, 
including the case where constant anomalous couplings are assumed for M wz < A FF , but 
invariant masses above App are ignored in deriving the sensitivity bounds |Kj . 

The bounds derived in this section are quite conservative. At the LHC, the limits can eas- 
ily be improved by 10 - 20% if W~Z + X production is included. Further improvements may 



result from using more powerful statistical tools than the simple \ test we performed 56 



IV. AMPLITUDE ZEROS, RAPIDITY CORRELATIONS, AND CROSS SECTION 

RATIOS 

Recently, it has been shown that the SM amplitude for qiq 2 — > W ± Z at the Born level 
exhibits an approximate zero at high energies, s 3> M§, located at || 

cos 6* = cos 6^ » ±- tan 2 9 W w ±0.1, (15) 
o 

where G* is the scattering angle of the Z boson relative to the quark direction in the WZ 
center of mass frame. The approximate zero is the combined result of an exact zero in the 
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dominant helicity amplitudes (±, =(=), and strong gauge cancellations in the remaining 
amplitudes. At high energies, only the (±, =f) and (0, 0) amplitudes remain non-zero in the 
SM. The existence of the zero in .M(±, =f) at cos©* ~ ±0.1 is a direct consequence of the 
contributing u- and t-channel fermion exchange diagrams and the left-handed coupling of 
the W boson to fermions. Unlike the W ±r y case with its massless photon kinematics, the 
zero has an energy dependence which is, however, rather weak for energies sufficiently above 
the W Z mass threshold. 

In this Section, we consider possible observable consequences of the approximate zero 
in WZ production in hadronic collisions and the impact of NLO QCD corrections on the 
relevant quantities. All numerical simulations are carried out using the parameters and cuts 
summarized in Sees. IIIA and IIIB. For the form factor, we again assume a dipole form factor 
(n = 2) with scale A FF = 1 TeV (see Eqs. (|7|) - @). Since the approximate amplitude zero 
in WZ production is similar in nature to the well-known radiation zero in W'j production, 
analogous strategies can be applied to search for observable signals. The radiation zero in 
W^f production leads to a pronounced dip in the rapidity distribution of the photon in the 
parton center of mass frame, da / dy*{^) pi|| . The approximate zero in the W Z amplitude is 



therefore expected to manifest itself as a dip in the corresponding y*{Z) distribution. Here, 

1 ^ f l+p z cos9* " 

2 11 1 1 -p z cos9* 



VIZ) = i in f i±£=£) , (16) 



and 



1 - 



AMI s 1 1/2 



(17) 



(s-M^ + Mf) 2 _ 

where, to lowest order, s = M^ z is the squared parton center of mass energy, and 9* the 
scattering angle of the Z boson with respect to the beam direction in the parton center of 
mass rest frame. For pp collisions the dip is centered at y*(Z) = 0. In pp collisions, the 
location of the minimum is determined by cos 0*> of Eq. ( |T5D , the average WZ invariant 
mass, and the fraction of events originating from sea quark collisions. As can be seen from 
Figs. § and ||, most of the cross section originates from the region \/§ = 200 - 250 GeV. 
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Valence quark collisions dominate at both, Tevatron and DiTevatron energies. The minimum 
of the y*(Z) distribution is therefore expected to occur at y*(Z) ~ ±0.06. 

The |t/*(Z)| distribution at the Tevatron and the LHC in Born approximation is shown in 



Fig. |16| The rapidity distribution of the Z boson in the parton center of mass frame at the 
DiTevatron is qualitatively very similar to that found at Tevatron energies and is therefore 
not shown. The SM |y*(Z)| distribution in the true WZ rest frame (dash-double-dotted 
curves) displays a pronounced dip at |y*(2')| = 0, which originates from the approximate 
amplitude zero. At Tevatron energies, the dip is quite significant. However, since the unob- 
servable longitudinal neutrino momentum can only be determined with a twofold ambiguity 
and, on an event to event basis, one does not know which solution is the correct one, both 
solutions have to be considered for each event. Assigning half of the event weight to each 
solution, the dip in the distribution using the reconstructed WZ rest frame is con- 

siderably filled in (solid lines). NLO QCD corrections further diminish the significance of 
the dip. 

In Fig. we also display the effect of non-standard WW Z couplings on the |y*(Z)| 
distribution (using the reconstructed WZ rest frame). The lower dotted, dashed, and dot- 
dashed lines display the \y*(Z)\ distribution for Ak° = +1 (+1), Ag° = +0.5 (+0.25), and 
A = +0.5 (+0.25), whereas the upper curves show the |t/*(Z)| spectrum for Aft = —1 
(-1), Ag® = -0.5 (-0.25), and A = -0.5 (-0.25) at the Tevatron (LHC). Non-standard 
WW Z couplings eliminate the approximate amplitude zero and, in general, tend to fill in 
the dip. However, due to the relatively strong interference between standard and anomalous 
contributions to the helicity amplitudes for An and Ag° at low energies, the dip may even 
become more pronounced for certain (positive) values of these two couplings at the Tevatron 
(see the lower dotted line in Fig. |16|a). 



From Fig. [16] it is obvious that the dip signaling the approximate amplitude zero in 
qiq 2 — > WZ will be difficult to observe in the |y*(Z)| distribution. In Wj production, 
correlations between the rapidities of the photon and the charged lepton originating from 



the W decay offer better access to the SM radiation zero than the y*(^y) distribution |L3 
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Knowledge of the neutrino longitudinal momentum, p L [y\ is not required in determining 
these correlations, and thus event reconstruction problems originating from the two possible 
solutions for Pjj^y) are automatically avoided. In 2 — > 2 reactions rapidity differences are 
invariant under boosts, Ay (7, W) = 2/(7)— y(W) = y*{y)—y*(W). One therefore expects the 
rapidity difference distribution, da/dAy(j, W), to exhibit a dip signaling the SM radiation 
zero. In W ±r y production, the dominant W helicity in the SM is Xw = ±1 ||57|| , implying 
that the charged lepton from W — > £v tends to be emitted in the direction of the parent 
W, and thus reflects most of its kinematic properties. As a result, the dip signaling the SM 
radiation zero manifests itself also in the Ay(j,£) = 2/(7) — y(£) distribution. 

The corresponding y(Z)—y(£\) distribution for W + Z production in the Born approxima- 
tion is shown in Fig. [L7| (solid line). Analogous to the situation encountered in qiq 2 — > Wj, 
the approximate zero in the WZ amplitude leads to a dip in the y(Z) — y(W) distribu- 
tion [llj], which is located at y(Z) — y(W) ~ ±0.12 (= 0) for W Z production in pp (pp) 



collisions. However, in contrast to W'-f production, none of the W helicities dominates in 
WZ production ]57| . The charged lepton, £1, originating from the W decay, W — > iiVi, 
thus only partly reflects the kinematical properties of the parent W boson. As a result, a 
significant part of the correlation present in the y(Z) — y{W) spectrum is lost, and only a 
slight dip survives in the SM y{Z) —y(£i) distribution. Due to the non-zero average rapidity 
difference between the lepton £\ and the parent W boson, the location of the minimum of 
the y(Z) — y{£\) distribution in pp collisions is slightly shifted to y(Z) — y(£i) ~ 0.5. The 
y(Z) —y(£i) distribution at the DiTevatron is qualitatively very similar to that obtained for 
pp collisions at y/s = 1.8 TeV (see Fig. [18]). 

The significance of the dip in the y(Z) —y(£i) distribution depends to some extent on the 
cut imposed on p T {£\) and the missing transverse momentum. Increasing (decreasing) the 
cut on p T {£\) (jf) T ) tends to increase the probability that £\ is emitted in the flight direction 
of the W boson, and thus enhances the significance of the dip. If the j) T > 50 GeV cut at 
the LHC could be reduced to 20 GeV, the dip signaling the approximate zero in the WZ 
production amplitude would be strengthened considerably. 
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Although the Z boson rapidity can readily be reconstructed from the four momenta of 
the lepton pair £^£2 originating from the Z decay, it would be experimentally easier to 
directly study the rapidity correlations between the charged leptons originating from the 
Z -> it £2 and W -> l x vx decays. The dotted lines in Fig. [T7| show the SM y(£j) - y(£f) 
distribution for W + Z production in the Born approximation. Because none of the Z or 
W helicities dominates []57j in qiq 2 — > WZ, the rapidities of the leptons from W and Z 
decays are almost completely uncorrelated, and essentially no trace of the dip signaling the 
approximate amplitude zero is left in the yil^) — y(£f) distribution. The y(l%) — y{£f) 
spectrum almost coincides with the y(l% ) — y{£f) distribution and is therefore not shown. 

In Figs. and [19] we show the influence of NLO QCD corrections and non-standard 
WW Z couplings (at NLO) on the Ay(Z, £1) = y(Z) —y(£i) spectrum. At Tevatron energies, 
the shape of the distribution is seen to be hardly influenced by the 0(as) QCD corrections. 
At the DiTevatron, the significance of the dip is slightly reduced. At LHC energies, the dip 
is completely eliminated by the inclusive QCD corrections. The NLO WZ + jet Ay(Z, £\) 
distribution, however, is very similar to the leading order rapidity difference distribution 



(see Fig. [T9|a). 

The effect of anomalous WW Z couplings on the NLO Ay(Z,£i) distribution is exem- 
plified by the dashed and dotted lines in Fig. [18] and in Fig. |l9|b. Similar to the situation 
encountered in the |t/*(Z)| distribution, the dip in the Ay(Z, £{) distribution at Tevatron and 
DiTevatron energies may be more pronounced than in the SM for certain (positive) values of 
Ak°. The shape of the Ay(Z, £\) distribution is seen to be quite sensitive to the sign of Ak° 
(Fig. 18|) . The same behaviour is observed for Ag®, whereas positive and negative values of 



A lead to a very similar Ay(Z,£i) distribution. In general, non-standard WW Z couplings 
tend to fill in the dip. In order not to overburden the figures, curves for A and Ag\ are not 
shown in Fig. [H| If deviations from the SM prediction were to be observed, it would be 
difficult to determine the sign of an anomalous coupling from the shape of the WZ invari- 
ant mass distribution, the cluster transverse mass spectrum, or the p T (Z) distribution. For 
Ak° and Ag®, the pronounced difference in shape of the Ay{Z 1 £\) distribution for positive 

30 



and negative values may aid in determining the sign. The influence of non-standard WWZ 
couplings on the exclusive NLO WZ + jet distribution is shown in Fig. |T9|b. Curves are 
only shown for positive values of the anomalous couplings. 



The error bars associated with the solid curves in Figs. [18] and |19| a indicate the expected 
statistical uncertainties for an integrated luminosity of 10 fb" 1 at the Tevatron and DiTe- 
vatron, and for J£dt = 100 fb _1 at the LHC It appears that the approximate zero in the 
WZ amplitude will be rather difficult to observe in the Ay(Z,£{) distribution. However, if 
the LHC is operated below its design luminosity of £ = 10 34 cm -2 s _1 , it may be possible to 
reduce the j> T cut. As mentioned above, the significance of the dip in the Ay(Z, t\) distribu- 
tion increases if the missing transverse momentum cut is lowered. The smaller collision rate 
is at least partially compensated by the larger total cross section for the reduced f T cut. It 
is thus possible that the conditions to detect the dip in the Ay{Z y £\) distribution improve 
if the LHC is operated below its design luminosity. However, more detailed simulations are 
required before definite conclusions can be drawn. 

As an alternative to rapidity correlations, cross section ratios can be studied. Many ex- 
perimental uncertainties, for example those associated with the lepton detection efficiencies, 
or the uncertainty in the integrated luminosity are expected to cancel, at least partially, 
in cross section ratios. In Ref. flTSl , the ratio of the Z7 to Wj cross sections, TZz-y/w-y, 
was shown to reflect the radiation zero present in the SM Wj helicity amplitudes. Due to 
the radiation zero, the W'j cross section is reduced in the central rapidity region. With 
increasing photon transverse momenta, events become more and more central in rapidity. 
The reduction of the cross section at small rapidities originating from the radiation zero 
thus becomes more pronounced at high Pxil)- This causes the photon transverse momentum 
distribution of q\(j2 — > to fall significantly faster than the p T {^) spectrum of qq — > Z7 
where no radiation zero is present. As a result, IZz-y/w-y increases rapidly with the minimum 
transverse momentum of the photon. 

In Fig. EO we study the cross section ratio 
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B\Z -» l+f ) <r(ZZ) _ B(Z -» gT) <r(ZZ) 

zz/wz B (Z ^£+£-)B(W ^£u)a(W ± Z) B(W ^ £u) a(W ± Z) ' 1 J 

as a function of the minimum transverse momentum of the Z boson, p™ in . To calculate 



the ZZ cross section, we use the results of Ref. |15| and assume the SM to be valid. The 
ZZ helicity amplitudes do not exhibit any zeros, whereas the SM WZ amplitude shows 
an approximate zero in the central rapidity region. The situation is thus qualitatively very 
similar to that encountered in the ratio of Z7 to cross sections, and one expects TZ Z z/wz 



to grow with p™ in . Figure |20|a demonstrates that, at Tevatron energies, TZzz/wz indeed rises 
quickly for p™ n > 100 GeV in the SM, indicating the presence of the approximate zero in the 
WZ amplitude. For smaller values of the minimum Z boson transverse momentum, TZzz/wz 
is approximately constant. In the low p™ in region, the shape of the p T (Z) distribution is 
dominated by Z mass effects which are similar in both processes. The cross section ratio at 
next-to-leading order differs only by about 10% from the LO ratio. 

At LHC energies, the situation is more complex. For p™ m < 100 GeV, TZzz/wz drops 
sharply due to the large f T cut imposed, which significantly suppresses the WZ cross section. 
While the cross section ratio slowly rises with p™ m for p™ m > 100 GeV at leading order, 
TZzz/wz continues to decrease if inclusive NLO QCD corrections are taken into account 
(Fig. pO[b). The relatively slower rise of TZzz/wz at LO at the LHC is due to the larger 
fraction of the cross sections originating from sea quark collisions, and the different x-ranges 
probed at the Tevatron and LHC. For p^ m (Z) = 1 TeV, the inclusive NLO cross section 
ratio is about a factor 3 smaller than TZzz/wz at leading order. At large values of the Z 
boson transverse momentum, the QCD corrections to WZ production at LHC energies are 



substantially larger than in the ZZ case |TJj , resulting in a large discrepancy between the LO 
and NLO prediction for TZzz/wz- in contrast to the situation encountered at the Tevatron, 
higher order QCD corrections completely blur the signal of the approximate amplitude zero 
in the WZ channel. Their size, however, can be substantially reduced by imposing a zero jet 
requirement (see Sec. HIE and Ref. JL5]]). The result for the ZZ + jet to W ± Z + jet cross 
section ratio at the LHC is given by the dotted line in Fig. pO[b. With a jet veto imposed, the 
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NLO ZZ to WZ cross section ratio rises with the minimum Z boson transverse momentum 
for p™ m > 100 GeV, and differs by at most 15% from the LO prediction. At the Tevatron, 
the NLO 0-jet cross section ratio virtually coincides with the ratio obtained at LO. 

The dot-dashed curve in Fig. |20], finally, shows the ZZ to WZ cross section ratio for 
Ak° = +1, illustrating the behaviour of IZzz/wz m presence of anomalous WW Z cou- 
plings. At the Tevatron (Fig. |20|a), the dot-dashed curve has been calculated taking into 
account inclusive 0(a s ) QCD corrections. At LHC energies (Fig. pOlb) the NLO ZZ to 
WZ cross section ratio is plotted with a jet veto included. Non-standard couplings lead 
to an enhancement of the WZ cross section, in particular at large values of Pp(Z) and, at 
the Tevatron, IZzz/wz decreases with p™ in . Due to the form factor parameters assumed 
(n = 2 and App = 1 TeV), the cross section ratio at the LHC displays a broad minimum 
at pp m (Z) « 300 GeV, and increases quickly at large values of p™ in . For larger values of 
App, and/or non-zero values of A g® or A , IZzz/wz rises more slowly, or may even decrease 
with p^ m (Z). In general, the ZZ to WZ cross section ratio as a function of the minimum 
Z boson transverse momentum differs substantially in shape from the SM prediction for 
IZzz/wz in presence of non-standard WW Z couplings. 

At the Tevatron, the limited number of ZZ and WZ events expected in the purely 
leptonic channels will unfortunately limit the usefulness of IZzz/wz- Even for an integrated 
luminosity of 10 fb _1 only a handful of events are expected for Pp(Z) > 150 GeV, and it will 
be very difficult to establish the growth with p™ m (Z) predicted by the SM. At the LHC, the 
statistical errors are expected to be much smaller, however, one can only hope to observe 
the rise of IZzz/wz signalling the presence of the approximate zero in the WZ channel if a 
0-jet requirement is imposed. Moreover, the rise of the cross section ratio is very slow, and 
for p™ m = 600 GeV only about 5 (2) purely leptonic W^Z (ZZ) events are expected for 
JCdt = 100 fb _1 . Combined, these effects will make it quite difficult to accurately determine 
the slope of IZzz/wz- 

A cross section ratio which suffers somewhat less from the small number of events ex- 
pected at the Tevatron, and which is less sensitive to QCD corrections at the LHC is the 
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ratio of WZ to W^y cross sections, 

Kwz/W " ~ B(W->ev)v(W±>y) -B(Z^££) (19) 

considered as a function of the minimum transverse momentum of the Z boson and photon, 
p™ m , respectively. TZwz/w-y measures the relative strength of the approximate zero in qiq 2 — » 

and the radiation zero in Wy production. Figure ^1] shows the ratio TZwz/w^ as a 
function of p™ m for the Tevatron (part a) and LHC (part b) center of mass energies. In 
obtaining TZwz/w^i we have considered both the electron and muon decay channels of the 
W and Z bosons. The W'y cross section in Fig. has been calculated using the results of 
Ref. ESI and the following cuts on the photon: 



p T (j) > 10 GeV, |r/( 7 )|<l, (Tevatron) 

p T ( 7 ) > 100 GeV, |r?( 7 )|<2.5, (LHC) (20) 

M T (£ T , i> T ) > 90 GeV, Ai?( 7 , £) > 0.7. 

At small values of p™ m , the WZ to cross section ratio rises very rapidly, due to 
the finite Z mass effects which dominate the shape of the Pt{Z) spectrum in this region 
for WZ production. For pf n > 100 GeV (200 GeV) at the Tevatron (LHC), K W z/ Wl is 
almost constant and independent of the center of mass energy, indicating that the radiation 
zero in qiq 2 — > W'y and the approximate amplitude zero in WZ production affect the 
respective photon and Z boson transverse momentum distribution in a very similar way. At 
the Tevatron, NLO QCD corrections reduce 'R-wz/w-y by about 10%. At LHC energies, the 
individual 0(a s ) QCD corrections are very large for both WZ and W'y production f49 



in 



particular at high transverse momenta (see Figs. |9] and [12]). In the cross section ratio, these 
large corrections cancel almost completely. For p™ m > 200 GeV, QCD corrections reduce 
TZwz/w-r by 20% or less. In contrast to the LO cross section ratio, which is completely flat 
for pf n > 200 GeV, K wz/Wl at NLO slowly rises with pf n at the LHC. 

Most theoretical models with non-standard WWZ couplings also predict anomalous 



WW'y couplings at the same time (see e.g., Ref. The effects of anomalous WWZ 
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and WW'-) couplings may cancel almost completely in TZwz/w^y if the WW Z and WW^f 
couplings are similar in magnitude and originate from operators of the same dimension. 
This is illustrated by the dot-dashed and dotted lines in Fig. ^T], which show TZwz/w-r at LO 
and NLO for Ak® = Ak° — — 1. Here the anomalous WW'j coupling Ak 7 is defined through 
an effective Lagrangian analogous to that of Eq. (|I|), and we assume equal form factor scales 
and powers (Aff = 1 TeV and n = 2) for An and Ak 7 . Both couplings correspond to 
operators of dimension four in the effective Lagrangian. Although the individual Pt{i) and 
Pt{Z) differential cross sections are enhanced by up to one order of magnitude (see e.g., 
Figs. [7| and |9|), IZwz/w-y agrees to better than 20% with the NLO SM cross section ratio for 
Ak° = Ak° = -I. 

At DiTevatron energies, the results for TZzz/wz and Hwz/w-y are qualitatively similar to 
those obtained for Tevatron and are therefore not shown. 

V. SUMMARY 

WZ production in hadronic collisions provides an opportunity to probe the structure 
of the WW Z vertex in a direct and essentially model independent way. Previous studies 
of pp — > W ± Z HUH have been based on leading order calculations. In this paper we have 
presented an 0(a s ) calculation of the reaction pp — > W ± Z+X — > ifv^i^+X for general, 
C and P conserving, WW Z couplings, using a combination of analytic and Monte Carlo 
integration techniques. The leptonic decays W — > t\i>\ and Z — > £2^2 have been included 
in the narrow width approximation in our calculation. Decay spin correlations are correctly 
taken into account in our approach, except in the finite virtual contribution. The finite 
virtual correction term contributes only at the few per cent level to the total NLO cross 
section, thus decay spin correlations can be safely ignored here. 

The Pt{Z) differential cross section is very sensitive to non-standard WW Z couplings. 
We found that QCD corrections significantly change the shape of this distribution at very 
high energies (see Fig. ^| and |T0]b). This shape change is due to a combination of destructive 
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interference in the WZ Born subprocess and a logarithmic enhancement factor in the qg 
and qg real emission subprocesses. The destructive interference suppresses the size of the 
WZ Born cross section and is also responsible for the approximate amplitude zero in qi<j2 — > 
WZ ||. The logarithmic enhancement factor originates in the high Pt(Z) [Pt(W)] region of 
phase space where the Z [W] boson is balanced by a high p T quark which radiates a soft W 
[Z] boson. The logarithmic enhancement factor and the large gluon density at high center 
of mass energies make the 0(a s ) corrections large for Pt{Z) 3> Mz- Since the Feynman 
diagrams responsible for the enhancement at large Pt{Z) do not involve the WW Z vertex, 
inclusive NLO QCD corrections to W^Z production tend to reduce the sensitivity to non- 
standard couplings. QCD corrections in WZ production thus exhibit the same features 
which characterize the 0(a s ) corrections in W^ production. 

At the Tevatron and DiTevatron, W Z production proceeds mainly via quark- ant iquark 
annihilation and, for the expected integrated luminosities (< 10 fb _1 ), large transverse 
momenta are not accessible. As a result, the sensitivity reduction in the high p T tail caused 
by the QCD corrections is balanced by the larger cross section at 0(a s ), and the limits 
derived from the NLO and LO Pt{Z) distribution are very similar (see Table I). At the 
LHC, however, where the qg luminosity is very high and the change in slope of the SM 
Pt{Z) distribution from QCD corrections is very pronounced, the sensitivity bounds which 
can be achieved are weakened by up to 40% (see Table II). 

The size of the QCD corrections at large p T {Z) may be reduced substantially, and a 
fraction of the sensitivity to anomalous WW Z couplings which was lost at the LHC may 
be regained, by imposing a jet veto, i.e., by considering the exclusive WZ + jet channel 
instead of inclusive WZ + X production. The improvement is equivalent to roughly a factor 
1.5 - 2.5 increase in integrated luminosity. The dependence of the NLO WZ + jet cross 
section on the factorization scale Q 2 is significantly reduced compared to that of the inclusive 
NLO WZ + X cross section. Uncertainties which originate from the variation of Q 2 will 
thus be smaller for sensitivity bounds obtained from the WZ + jet channel than for those 
derived from the inclusive NLO WZ + X cross section. At the Tevatron (DiTevatron), NLO 

36 



QCD corrections do not influence the sensitivity limits in a significant way. Nevertheless, it 
will be important to take these corrections into account when extracting information on the 
structure of the WWZ vertex, in order to reduce systematic and theoretical errors. 

At the Tevatron (DiTevatron) with fCdt = 10 fb _1 , taking into account all correlations 
between the different WWZ couplings, Ak° can be measured with 70 - 100% (50 - 60%) 
accuracy in WZ production in the purely leptonic channels, whereas the two other couplings 
can be determined with an uncertainty of 0.1 - 0.25. At the LHC with J Cdt = 100 fb , An 
can be determined with an uncertainty of about 10%, whereas Ag® and A can be measured 
to better than 0.01, with details depending on the form factor scale assumed (see Table II). 

The bounds listed in Tables I and II should be compared with the limits which can be 
obtained in other channels, and at e + e~ colliders. Assuming SM WW'y couplings, Ag® and 
A can be measured in pp — > W + W~ , W^Z — > i .ujj with a precision similar to that which 
can be achieved in the W ± Z — > ifvii^l^ mode, both at the Tevatron with 1 fb -1 and the 
TeV*. The limits which can be obtained for Ak° from the i^vjj final state are about a 
factor 4 better than those from double leptonic WZ decays The bounds which can be 
achieved for Ak°, Ag\ and A in e + e~ — > W + W~ at LEP II depend quite sensitively on the 
center of mass energy. For = 176 GeV and J Cdt = 500 pb _1 , the WWZ couplings can 
be measured with a precision of about ±0.5, if correlations between the three couplings are 



taken into account [53,58|. At a linear e + e collider with a center of mass energy of 500 GeV 



or higher, they can be determined with an accuracy of better than 0.01 [p3| , [59 . 

We also studied possible experimental signals of the approximate zero in the SM WZ 
amplitude. Unlike the situation encountered in W^ production where the radiation zero 
leads to a pronounced minimum in the photon-lepton rapidity difference distribution, the 
approximate amplitude zero in WZ production causes a slight dip only in the corresponding 
Ay(Z,£i) = y(Z) — y(£i) distribution. In W 7 production, the dominant W helicity is 
Xw = ±1, implying that the charged lepton from the decaying W boson tends to be emitted 
in the direction of the parent W boson, and thus reflects most of its kinematic properties. 
In contrast, none of the W helicities dominates in WZ production. The charged lepton 
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originating from the W boson decay, W — > t\V\, thus only partly reflects the kinematic 
properties of the parent W boson, which reduces the significance of the dip. At Tevatron 
and DiTevatron energies, higher order QCD corrections only negligibly influence the shape 
of the Ay(Z,£ 1 ) distribution. At the LHC, however, NLO QCD effects completely obscure 
the dip, unless a 0-jet requirement is imposed. 

Alternatively, cross section ratios can be used to search for experimental footprints of 
the approximate amplitude zero. We found that the ratio of ZZ to WZ cross sections as 
a function of the minimum Z boson transverse momentum, p™ m , increases with p™ in for 
values larger than 100 GeV. The increase of the ratio of ZZ to WZ cross sections is a 
direct consequence of the approximate zero. The ratio of WZ to W7 cross sections, on 
the other hand, is almost independent of the minimum p T of the Z boson and photon for 
sufficiently large values of p™ in , indicating that the approximate zero in WZ production and 
the radiation zero in Wj production affect the Z boson and photon transverse momentum 
distributions in a very similar way. QCD corrections have a significant impact on the ZZ to 
WZ cross section ratio at the LHC unless a jet veto is imposed, whereas they largely cancel 
in the WZ to W^ cross section ratio. 

Together with the Ay(Z,£ 1 ) distribution, the ZZ to W ± Z and WZ to W7 cross section 
ratios are useful tools in searching for the approximate amplitude zero in WZ production. 
However, for the integrated luminosities envisioned, it will not be easy to conclusively es- 
tablish the approximate amplitude zero in WZ production at the Tevatron or the LHC. 
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TABLES 

TABLE I. Sensitivities achievable at the 95% confidence level (CL) for the anomalous WW Z 
couplings A<7i, Ak°, and A in pp —* W ± Z + X -> ifvii^iz + X at a) leading order and b) 
next-to-leading order for the Tevatron, the TeV* (y/s = 1.8 TeV in both cases), and the DiTe- 
vatron (^/s = 3.5 TeV). The limits for each coupling apply for arbitrary values of the two other 
couplings. For the form factors we use Eqs. ([?]), (||), and (|9|) with n = 2 and Aff = 1 TeV. The 
cuts summarized in Sec. IIIB are imposed. 



coupling 



Tevatron 
j'Cdt = 1 ft)- 1 



a) leading order 

TeV* 
j'Cdt = 10 fb- 1 



DiTevatron 
fCdt = 10 fb" 1 



A 5 ? 



+0.52 
-0.29 

+1.9 
-1.4 

+0.34 
-0.37 



+0.26 
-0.10 

+1.0 
-0.7 

+0.15 
-0.15 



+0.15 
-0.05 

+0.5 
-0.5 

+0.08 
-0.07 



coupling 



b) next-to-leading order 
Tevatron TeV* 
jCdt = 1 fb" 1 jCdt = 10 fb" 1 



DiTevatron 
JCdt = 10 fb- 1 



A 5 ? 



Ak 



A 



+0.50 
-0.27 

+1.8 
-1.2 

+0.33 
-0.34 



+0.25 
-0.10 

+1.0 
-0.7 

+0.15 
-0.14 



+0.15 
-0.05 

+0.6 
-0.5 

+0.08 
-0.08 
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TABLE II. Sensitivities achievable at the 95% confidence level (CL) for the anomalous WWZ 
couplings A^, An , and A in pp -> W + Z + X -»■ # + X at the LHC (^i = 14 TeV). The 

limits for each coupling apply for arbitrary values of the two other couplings. For the form factors 
we use Eqs. (^), @, and (||) with n = 2. The cuts summarized in Sec. IIIB are imposed. In the 
NLO 0-jet case we have used the jet definition of Eq. (|i~4|). 





a) / Cdt = 


= 10 ft) -1 , 




= 1 TeV 




coupling 


Born appr. 






incl. NLO 


NLO 0-jet 




+0.108 






+0.127 


+0.115 


-0.036 






-0.044 


-0.038 


Ak° 


+0.53 
-0.46 






+0.62 
-0.59 


+0.60 
-0.49 


X° 


+0.058 






+0.065 


+0.063 


-0.057 






-0.073 


-0.062 




b) J Cdt = 


100 ftr 1 


, Aff 


= 1 TeV 




coupling 


Born appr. 






incl. NLO 


NLO 0-jet 




+0.082 






+0.096 


+0.085 


-0.014 






-0.018 


-0.016 


Ak° 


+0.17 
-0.34 






+0.24 
-0.41 


+0.20 
-0.35 


X° 


+0.038 






+0.042 


+0.036 


-0.036 






-0.048 


-0.038 




c) / Cdt = 


100 fbr 1 . 


Aff 


= 3 TeV 




coupling 


Born appr. 






incl. NLO 


NLO 0-jet 


A 5 ? 


+0.0164 






+0.0200 


+0.0188 


-0.0048 






-0.0066 


-0.0050 


Ak° 


+0.092 
-0.120 






+0.128 
-0.160 


+0.108 
-0.132 




+0.0084 






+0.0102 


+0.0092 


-0.0082 






-0.0100 


-0.0090 
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FIGURES 

FIG. 1. Feynman rule for the general WW Z vertex. The factor gwwz = ecot #vv is the WW Z 
coupling strength and Qw is the electric charge of the W boson. The vertex function Tp^ik, ki, 
is given in Eq. @. 

FIG. 2. The inclusive differential cross section for the reconstructed W Z mass in the reaction 
pp — > W + Z + X — » l\v\IQ £2 + X at -y/s = 1.8 TeV; a) in the Born approximation and b) includ- 
ing NLO QCD corrections. The curves are for the SM (solid lines), A = —0.5 (dashed lines), 
Ak° = —1.0 (dotted lines), and Ag® = —0.5 (dot-dashed lines). The cuts imposed are summarized 
in Sec. IIIB. 

FIG. 3. The inclusive differential cross section for the reconstructed WZ mass in the reaction 
pp — > W + Z + X — > vii.^ £2 + X at yfs = 3.5 TeV; a) in the Born approximation and b) includ- 
ing NLO QCD corrections. The curves are for the SM (solid lines), A = —0.5 (dashed lines), 
AkP = —1.0 (dotted lines), and Ag® = —0.5 (dot-dashed lines). The cuts imposed are summarized 
in Sec. IIIB. 

FIG. 4. The inclusive differential cross section for the reconstructed WZ mass in the reaction 
pp — ► W + Z + X — > l\ v\l~2 £2 + X at y/s = 14 TeV; a) in the Born approximation and b) including 
NLO QCD corrections. The curves are for the SM (solid lines), A = —0.25 (dashed lines), 
Ak° = —1.0 (dotted lines), and Ag® = —0.25 (dot-dashed lines). The cuts imposed are summarized 
in Sec. IIIB. 

FIG. 5. The inclusive NLO differential cross section for the reconstructed WZ mass 
in the reaction a) pp -> W + Z + X -> ifu^i^ + x at y/s = 1.8 TeV and b) pp -> W + Z + X 
— ► €[v\£~2 £2 + X at = 14 TeV. The curves are for the SM with reconstructed invariant mass 
(solid lines) and non-standard WW Z couplings (dotted, dashed, and dot-dashed curves) as listed 
on the figure. The lower (upper) lines apply for positive (negative) anomalous couplings. The 
dash-double-dotted line shows the true SM W Z invariant mass distribution. The cuts imposed are 
summarized in Sec. IIIB. 
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FIG. 6. The inclusive NLO differential cross section for the cluster transverse mass for a) 
pp -► W + Z + X -> Ifvilpz + x at Vs = 1-8 TeV and b) pp -► VF+Z + X £+i^+£T + X at 

= 14 TeV. The curves are for the SM (solid lines) and non-standard WW Z couplings (dotted, 
dashed, and dot-dashed curves) as listed on the figure. The cuts imposed are summarized in 
Sec. IIIB. 

FIG. 7. The inclusive differential cross section for the transverse momentum of the Z boson in 
the reaction pp -► W + Z + X -► l+v^q + X at v/s = 1.8 TeV; a) in the Born approximation and 
b) including NLO QCD corrections. The curves are for the SM (solid lines), A = —0.5 (dashed 
lines), Ak° = —1.0 (dotted lines), and Ag® = —0.5 (dot-dashed lines). The cuts imposed are 
summarized in Sec. IIIB. 

FIG. 8. The inclusive differential cross section for the transverse momentum of the Z boson in 
the reaction pp — ► W + Z + X — > l^vxl^l^ + X at \fs = 3.5 TeV; a) in the Born approximation and 
b) including NLO QCD corrections. The curves are for the SM (solid lines), A = -0.5 (dashed 
lines), Ak° = —1.0 (dotted lines), and Ag® = —0.5 (dot-dashed lines). The cuts imposed are 
summarized in Sec. IIIB. 

FIG. 9. The inclusive differential cross section for the transverse momentum of the Z boson in 
the reaction pp — > W + Z + X —> ifvii^^ + X at y/s = 14 TeV; a) in the Born approximation and 
b) including NLO QCD corrections. The curves are for the SM (solid lines), A = —0.25 (dashed 
lines), Ak° = —1.0 (dotted lines), and Ag® = —0.25 (dot-dashed lines). The cuts imposed are 
summarized in Sec. IIIB. 

FIG. 10. Ratio of the NLO to LO differential cross sections of the Z boson trans- 
verse momentum as a function of p T {Z) for a) pp — > W + Z + X — > tyVxt^l^ + X and b) 
pp W + Z + X -> lXv x l^i^ + X at = 14 TeV. In part a) the solid (dotted) and dashed 
(dot-dashed) lines give the ratio in the SM and for A = —0.5 at the Tevatron (DiTevatron), 
respectively. In part b) the solid and dashed curves show the ratio for the SM and A = —0.25 at 
the LHC. The cuts imposed are summarized in Sec. IIIB. 
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FIG. 11. The differential cross section for the Z boson transverse momentum in the reaction 
pp -> W + Z + X -> ifvii^q + X at y/a = 1.8 TeV in the SM. a) The inclusive NLO differential 
cross section (solid line) is shown, together with the 0(a s ) 0-jet (dotted line) and the (LO) 1-jet (dot 
dashed line) exclusive differential cross sections, using the jet definition in Eq. (||). b) The NLO 
WZ + jet exclusive differential cross section (dotted line) is compared with the Born differential 
cross section (dashed line). The cuts imposed are summarized in Sec. IIIB. 

FIG. 12. The differential cross section for the Z boson transverse momentum in the reaction 
pp -> W + Z + X -> i\v x i^i 2 + X at = 14 TeV in the SM. a) The inclusive NLO differential 
cross section (solid line) is shown, together with the 0(a s ) 0-jet (dotted line) and the (LO) 1-jet (dot 
dashed line) exclusive differential cross sections, using the jet definition in Eq. (0). b) The NLO 
WZ + jet exclusive differential cross section (dotted line) is compared with the Born differential 
cross section (dashed line). The cuts imposed are summarized in Sec. IIIB. 

FIG. 13. The total cross section for pp — » W + Z + X — > ifuxi^t^ + X in the SM versus the 
scale Q; a) at the Tevatron and b) at the LHC. The curves represent the inclusive NLO (solid lines), 
the Born (dot-dashed lines), the LO 1-jet exclusive (dashed lines), and the NLO 0-jet exclusive 
(dotted lines) cross sections. The cuts imposed are summarized in Sec. IIIB. For the jet definitions, 
we have used Eqs. (ll~3| ) and (|l4|). 

FIG. 14. Limit contours at the 95% CL for pp -> W ± Z + X -> tfv^t^ + X derived from the 
inclusive NLO Pt(Z) distribution. Contours are shown in three planes: a) the AkP - A plane, b) 
the AkP - Ag® plane, and c) the Ag® - X° plane. The solid and dashed lines give the results for 
the Tevatron (y/s = 1.8 TeV) with fCdt = 1 ftr 1 and fCdt = 10 fb" 1 , respectively. The dotted 
curve shows the result obtained for the DiTevatron (y/s = 3.5 TeV) with fCdt = 10 fb . The 
cuts imposed are summarized in Sec. IIIB. 
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FIG. 15. Limit contours at the 95% CL for pp -► W + Z + X -> i+u^+q + X at y/s = 14 TeV 
derived from the p T {Z) distribution. Contours are shown in three planes: a) the An - X° plane, 
b) the Ak°- Agi plane, and c) the Ag®- \° plane. The solid and dashed lines give the inclusive 
NLO and LO results, respectively, for / Cdt = 10 fb -1 . The dotted and dot-dashed curves show the 
results obtained from the exclusive NLO WZ + jet channel for integrated luminosities of 10 fb -1 
and 100 fb -1 , respectively. The cuts imposed are summarized in Sec. IIIB. 

FIG. 16. Rapidity spectrum of the Z boson in the WZ rest frame in the Born approximation for 
a) pp -► W + Z + X -> l+v^+q + X at <Ji = 1.8 TeV and b) pp -► W + Z + X -> l\v x l^q + X 
at = 14 TeV. The curves are for the SM with reconstructed WZ rest frame (solid lines) and 
non-standard WW Z couplings (dotted, dashed, and dot-dashed curves) as listed on the figure. The 
lower (upper) lines apply for positive (negative) anomalous couplings. The dash-double-dotted line 
shows the true SM |y*(Z)| distribution. The cuts imposed are summarized in Sec. IIIB. 

FIG. 17. SM rapidity difference distributions in the Born approximation for a) pp — > 
W + Z + X ^ l+v^+q +X at v/i = 1.8 TeV and b) pp -► W + Z + X -> i+v^+q + X at 
y^s = 14 TeV. The curves are for y(Z) — y((f) (solid lines) and yil^) — y{i\) (dotted lines). 
The cuts imposed are summarized in Sec. IIIB. 

FIG. 18. The differential cross section for the rapidity difference Ay(Z,£i) for pp — > 
W + Z + X -> t+vxl+q, + X a) at ^ = 1.8 TeV and b) at y/s = 3.5 TeV. The solid and dot-dashed 
curves show the inclusive NLO and the LO SM prediction, respectively. The dashed and dotted 
lines give the results for AkP = +1 and Ak° = — 1, respectively. The error bars associated with the 
solid curves indicate the expected statistical uncertainties for an integrated luminosity of 10 fb -1 . 
The cuts imposed are summarized in Sec. IIIB. 
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FIG. 19. The differential cross section for the rapidity difference Ay(Z,£i) at y/s = 14 TeV for 
a) pp -> W + Z + X -> l\v x l\q + X in the SM and b) for pp -» VF+Z + jet -> #*i£^ + jet 
at NLO. In part a) the dotted and dashed curves show the inclusive NLO and the LO SM prediction, 
respectively, while the solid curve gives the prediction for the SM NLO W Z + jet case. The error 
bars associated with the solid curve indicate the expected statistical uncertainties for an integrated 
luminosity of 100 fb _1 . In part b) the curves are for the SM (solid line), Ak° = +1 (dotted 
line), A = +0.25 (dashed curve), and Ag\ = +0.25 (dot-dashed curve). The cuts imposed are 



summarized in Sec. IIIB. For the jet definition, we have used Eq. (14). 



FIG. 20. The ratio TZ ZZ/WZ = B(Z t+l~) a(ZZ)/B(W -► £v)a(W ± Z), I = e, /j, as a 
function of the minimum transverse momentum of the Z boson, p™ m , at a) the Tevatron and 
b) the LHC. The solid and dashed line show the inclusive NLO and the LO result for the SM, 
respectively. The dotted line in b) gives the SM cross section ratio at NLO if a 0-jet requirement is 
imposed. The dot-dashed line displays IZzz/wz f° r = +1- In part a) this curve is calculated 
taking into account inclusive NLO QCD corrections, whereas in part b) the dot-dashed curve is 
for the NLO 0-jet cross section ratio. The cuts imposed are summarized in Sec. IIIB. For the jet 



definition, we have used Eq. (14) 



FIG. 21. The ratio TZ WZ/W ^ = B{Z -> t+tr) a(W ± Z)/a(W ± j), £ = e, fi, as a function of the 
minimum transverse momentum of the Z boson and photon, p™ in , respectively, at a) the Tevatron 
and b) the LHC. The solid and dashed lines show the inclusive NLO and the LO result for the 
SM, respectively. The dotted and dot-dashed lines display the inclusive NLO and LO WZ to Wj 
cross section ratio for Ak!J = AkP = — 1. Here, Ak 7 is the anomalous WWj coupling defined in 
an analogous way to An [see Eq. (|l])]. The cuts imposed are summarized in Sec. IIIB and Eq. (| 
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